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Abstract 

<N 

^ The theory of dissipativity has been primarily developed for controllable systems/behaviors. For 

various reasons, in the context of uncontrollable systems/behaviors, a more appropriate definition of 
dissipativity is in terms of the dissipation inequality, namely the existence of a storage function. A 
storage function is a function such that along every system trajectory, the rate of increase of the storage 
function does not exceed the power supplied. While the power supplied is always expressed in terms of 
only the external variables, whether or not the storage function should be allowed to depend on only the 
external variables and their derivatives or also unobservable/hidden variables has various consequences 
on the notion of dissipativity: this paper thoroughly investigates the key aspects of both cases, and 
also proposes another intuitive definition of dissipativity. 

Wc first assume that the storage function can be expressed in terms of the external variables 
I— ' and their derivatives only and prove our first main result that, assuming the uncontrollable poles are 

unmixed, i.e. no pair of uncontrollable poles add to zero, and assuming a strictness of dissipativity 
at the infinity frequency, the dissipativities of a system and its controllable part are equivalent; in 
other words once the autonomous subsystem satisfies a Lyapunov equation solvability-like condition, 
it does not interfere with the dissipativity of the system. We also show that the storage function in 
this case is a static state ftmction. This main result proof involves new results about solvability of 
the Algebraic Riccati Equation, and uses techniques from Indefinite Linear Algebra and Hamiltonian 
^-H matrix properties. 

We then investigate the utility of unobservable/hidden variables in the definition of storage function: 
. ^ we prove that lossless uncontrollable behaviors are ones which require storage function to be expressed 

in terms of variables that are unobservablc from the external variables. 
^ Wc next propose another intuitive definition: a behavior is called dissipative if it can be embedded 

in a controllable dissipative super-behavior. We show that this definition imposes a constraint on 
the number of inputs and thus explains unintuitive examples from the literature in the context of 
lossless/orthogonal behaviors. These results are finally related to RLC realizability of passive networks, 
specifically to the nonrealizability of the nuUator one-port circuit using RLC elements. 



<N 
O 



1 Introduction 

The theory of dissipativity for linear dynamical systems helps in the analysis and design of control systems 
for several control problems, for example, LQR/LQG control, "K^o, synthesis of passive systems, and 
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optimal estimation problems. When dealing with LTI systems, it is straightforward to define dissipativity 
for controllable systems due to a certain property of such systems that their compactly supported system 
trajectories are, loosely speaking, 'dense' in the set of all allowed trajectories. However, this is not the 
case for uncontrollable systems, and this situation is the central focus of this paper. We elaborate more 
on this point when we define dissipativity and review equivalent conditions for controllable systems in 
Section H 

In this paper, we use a less-often-used definition of dissipativity for systems, possibly uncontrollable, 
and generalize key results using some techniques from indefinite linear algebra (see [10]) for solving 
Algebraic Riccati Inequalities in the context of an uncontrollable state space system. Like in [TJ [12] , we 
define a system as dissipative if there exists a storage function that satisfies the dissipation inequality 
for all system trajectories. The existential aspect of this definition raises key issues that this paper deals 
with. 

The main result we show is that if the uncontrollable poles of an LTI system are such that no two 
of them add to zero, and if the controllable subsystem strictly dissipates energy at frequency equal to 
infinity, then the dissipativity of the controllable subsystem is equivalent to the system's dissipativity. 
We also show that, using the concatenability axiom of the state, the energy stored in a system is a static 
function of the state variables. Further, we also show that this state is 'observable' from the external 
variables, i.e. the state is a linear combination of the external variables and possibly their derivatives. 
This is intuitively expected in view of the fact that energy exchange between the system and its ambience 
takes place through the external variables. However, it appears that this may not be the case for lossless 
systems, i.e. systems that don't dissipate any energy, nor contain a source within. 



2 Preliminaries 

In this section we include various definitions about the behavioral framework for studying dynamical 



systems (Subsection 2.1) and then introduce background results about dissipative systems (Subsection 
2.2). Subsection |2.3| contains brief notation about indefinite linear algebra from [10]. For this paper, 
M denotes the set of all real numbers and M[^] the set of polynomials in the indeterminate ^ and real 
coefficients; matrices and polynomial matrices are denoted the obvious way. We use • to leave a row 
dimension unspecified, for example, ]R*^'^[,^]. The space of infinitely often differentiable functions from M 
to say M" is denoted by C°°(M, M") and D denotes the set of compactly supported functions within this 
space. 



2.1 The behavioral approach 

When dealing with linear differential systems, it is convenient to use polynomial matrices for describing 
a differential equation. Suppose Rq, Ri . . . Rjy are constant matrices of the same size such that 

Rqw + Ri^w + R2^w . . . Rn^w = 

is a linear constant coefficient ordinary differential equation in the variable w. We define the polynomial 
matrix R{^) := Rq + RiS, + ■ ■ ■ Rn^^ , and represent the above differential equation as i?(4-)u) = 0. 
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A linear differential behavior, denoted by 5S, is defined as the set of all infinitely often differentiable 
trajectories that satisfy a system of ordinary linear differential equations with constant coefficients, i.e.. 



where R{(,) is a polynomial matrix having w nmnber of columns, i.e., R{(,) G M*^"'[^]. We denote the set 
of all such linear differential behaviors with w number of variables by The linear differential behavior 
!B G given by the above definition can also be represented by 53 = keri?(^). One views the rows of 
R as differential equations that the variable w has to satisfy in order for a trajectory w{t) to be in the 
behavior 55. The matrix R is not unique and one can use elementary row operations to modify R and 
this does not change the set of solutions 53: this thus allows assuming without loss of generality that R 
has full row rank (see [2¥j). The number of inputs of !B is defined as w — rank (R) and is called the input 
cardinality. This integer depends only on *B and not on the R used to define it; m(5S) denotes the number 
of inputs. 

An important fundamental concept is controllability of a system. A behavior 53 = kei R{^)w is said 
to be controllable, if for every wi and W2 £ there exist t^a G 53 and r > such that 



The set of all controllable behaviors with w variables is denoted as £"cont- This patchability definition of 
controllability is known to have the traditional Kalman state-space definition of controllability as a special 
case in |24| . Further, it is shown that 53 = kei R{^)w is controllable if and only if R{X) has constant 
rank for all A G C. It is also shown in |24j that 53 is controllable if and only if it can be defined as 



where M(^) G M"^'"[^]. This representation of 53 G 2,^ont known as an image representation. The 
variable £ is called a latent variable: these are auxiliary variables used to describe the behavior; we 
distinguish the variable w as the manifest variable, the variable of interest. It is known that 53 G ^^ont 
always allows an image representation with M(^) such that M(A) has full column rank for every A G C. 
This kind of image representation is known as an observable image representation. In this paper, unless 
otherwise stated explicitly, we assume the image representations are observable. The use of the term 
'observable' is motivated by the fact that the variable i is observable from the variable w. This notion is 
defined as follows. 

For a behavior !B with variables w and £, we say i is observable from w if whenever (w, ii) and (w, I2) 
both are in *B, we have ii = £2- Observability of i from w in a behavior 53 implies that there exists a 
polynomial matrix F(^) such that £ = F{^)w for all w and i in the behavior. 

We now define relevant notions in the context of uncontrollable behaviors. For a behavior 53, possibly 
uncontrollable, the largest controllable behavior contained in 53 is called the controllable part of 53, and 
denoted by ^cont- An important fact about the controllable part of 53 is that m(53cont) = ni(53). The set 
of complex numbers A for which R{\) loses rank is called the set of uncontrollable modes and is denoted 
by Aun- For a detailed exposition on behaviors, controllability and observability we refer the reader to 



53 := G (J:°°(M,M") | R{—)w = 0} 




wi (t) for all t ^ 0, 
W2{t) for all t ^ r. 



53 := {u; G e:°°(M,M") | there exists an^ G (j:°°(M,M'") 
such that ID = M(^)£}, 



3 



2.2 Quadratic Differential Forms and dissipativity 

The concept of Quadratic Differential Forms (QDF) (see |20j) is central to this paper. Consider a two 
variable polynomial matrix with real coefficients, rj) := ^jkC^rj^ € M"^" [C, rj], where S M"^", 
The QDF Q$ induced by ^> is a map : e:°°(M,]R") — > (t°°{R,R) defined as 

When dealing with quadratic forms in w and its derivatives, we can assume without loss of generality 
that = C)- such a <I> is called a symmetric two variable polynomial matrix. A quadratic 

form induced by a real symmetric constant matrix S £ M"^" is a special QDF. We frequently need the 
number of positive and negative eigenvalues of a nonsingular, symmetric matrix S: they are denoted by 
cr+(S') and cr_(5) respectively. 

For a two variable polynomial matrix ^{(^,r]) € R"^" [Ci^]i define the single variable polynomial 
matrix 9$ by d^{S,) = ^{-^,0- 

Consider T, € M'"^"', a symmetric nonsingular matrix. A behavior 5S G is said to be dissipative 
with respect to the supply rate S (or S-dissipative) if there exists a QDF such that 

itQ^iw) ^ w'^^w for allwG'B . (1) 

The QDF Qi^, is called a storage function: it signifies the energy stored in the system at any time 
instant. The above inequality is called the dissipation inequality. A behavior !B is called S-lossless if the 
above inequality is satisfied with an equality for some QDF Qijf. Notice that the storage function plays 
the same role as that of Lyapunov functions in the context of autonomous systems; the notion of storage 
functions is a generalization to non-autonomous systems of Lyapunov functions, as pointed in |20j . The 
following theorem from [20] applies to controllable behaviors. 

Proposition 2.1 Consider 55 G -C"cont O'^d let w = M{-j^)£ be an observable image representation. 
Suppose T, G M^^"^ is symmetric and nonsingular. Then, the following are equivalent. 

1. There exists a QDF Qq, such that inequality ([T]) is satisfied for all w 

2. f^w'^Tiwdt ^ for all w £ ^ (I'D, the compactly supported trajectories in 53. 

3. M^{-juj)T.M{ju}) ^ for all w G M. 

The significance of the above theorem is that, for controllable systems, it is possible to verify dissi- 
pativity by checking non-negativity of the above integral over all compactly supported trajectories: the 
compact support signifying that we calculate the 'net power' transferred when the system starts 'from 
rest' and 'ends at rest'. The starting and ending 'at rest' ensures that for linear systems there is no 
internal energy at this time. The absence of internal energy allows ruling out the storage function from 
this condition: in fact, this is used as the definition of dissipativity for controllable systems. The same 
cannot be done for uncontrollable systems due to the compactly supported trajectories not being dense in 
the behavior (see pL3]). An extreme case is an autonomous behavior, i.e. a behavior which has m(?B) = 0: 
while the zero trajectory is the only compactly supported trajectory, the behavior consists of exponentials 
corresponding to the uncontrollable poles of 'B. The issue of existence of storage functions is elaborated 
in [20l Remark 5.9] and in text following [12l Proposition 3.3]. 
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2.3 States 



The state variable is defined as a latent variable that satisfies the property of state, that is, if {wi,xi), 
{w2,X2) G ^/ and a;i(0) = X2(0), then the new trajectory {w,x) formed by concatenating (wi,a:i) and 
{w2, X2) at t = 0, i.e., 

{wi,xi){t) for all t^O 

[W, X){t) = 

{w2,X2){t) for all t > 0, 

also satisfies the system equation of *Bj in a distributional sense |21j. It is intuitively expected that a 
variable x has the state property if and only if w and x satisfy an equation that is at most first order 
in X and zeroth order in w: see [15] for precise statement formulation and proof. When w is partitioned 
into w = {wi,W2), with wi as the input and W2 as the output, then 03/ admits the more familiar 
input/state/output (i/s/o) representation as 

—x = Ax + B'Wi, W2 = Cx + Dwi. (2) 

One can ensure that x is observable from w; this is equivalent to conventional observability of the pair 
(C, A). While such a state space representation is admitted by any 5S, the pair [A, B) is controllable (in 
the state space sense) if and only if !B is controllable (in the behavioral sense defined above). 

2.4 Indefinite linear algebra 

In this paper, we use certain properties of matrices that are self-adjoint with respect to an indefinite inner 
product. We briefly review self-adjoint matrices and neutral subspaces (see jlQj). Let P G c^x" an 
invertible Hermitian matrix. This defines an indefinite inner product on C" by {Px, x) := x*Px, where x* 
is the complex conjugate transpose of the vector x £ C". For a complex matrix A, the complex conjugate 
transpose of A is denoted by A*. 

Consider matrices A and P € C"^" with P invertible and Hermitian. The P-adjoint of the matrix A, 
denoted by A^*\ is defined as A^*^ = P~^A*P. The matrix A is said to be P-self-adjoint if A = A^*\ i.e., 
A = P-^A*P. 

Since P is not sign-definite in general, the sign of (Px, x) is zero, positive or negative depending on 
the vector x. A subspace M C is said to be P-neutral if the inner product {Px, x) = for all x G M. 

3 Dissipativity of uncontrollable behaviors 

In this paper, we deal with systems which satisfy the dissipativity property, i.e. net energy is directed 
inwards along every system trajectory. As elaborated below, the 'total' aspect of the energy involves 
an integral, thus bringing in the initial and final conditions of the trajectory being integrated. The 
convenience of starting-from-rest and ending-at-rest applies to only controllable systems as we will review 
soon. The notion of storage function helps in formulating the dissipation property as an inequality to be 
satisfied at each time-instant. A central issue in this paper is what variables should the storage function 
be allowed to depend on. We explore dependencies on a latent variable i, or on a state variable x, or on 
the manifest variable w. this is indicated in the storage function subscript; the following definition has 
appeared in several works, see [201 Ell [2lj , for example. 
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Definition 3.1 Let S G M'^^"' he a nonsingular symmetric matrix, inducing the supply rate w'^TjW. 
Consider a behavior 5S G £" with manifest variables w and latent variable £, with the corresponding full 
behavior 5Sf^jj. For the behavior let x be a state variable with the corresponding full behavior 55 j^^. 
With respect to the supply rate T,, the behavior OS is said to be dissipative if there exists a quadratic 
differential form such that 

^tQ^M) ^ w^^w for all {w,£)€ 'B^. (3) 

1. The function Q^^, a quadratic function of d. and its derivatives, is called a storage function. 

2. A storage function is said to be an observable storage function if the latent variable I is ob- 
servable from the manifest variable w. In this case, there exists a storage function Q^^ such that 
^Q^^iw) ^ vu'^Y.uj for allw^'^B. 

3. A storage function Qq,^ (x) is said to be a state function if Q<i,^ (x) is equal to x'^Kx for some 
constant matrix K . 

In this paper, we study dissipativity with respect to a constant, nonsingular, symmetric S G M"^". 
It is known (see Remark 5.11] and j21l Proposition 2]) that S-dissipativity of a behavior !B implies 
that the input cardinality of !B cannot exceed (T+(S), i.e., ni(5S) ^ cr+(S). In this context it is helpful 
to perform a coordinate transformation in the variables w so that S is a diagonal matrix consisting of 
only +l's and — I's along the diagonal. Moreover, there exists an input/output partition such that all the 
inputs correspond to +l's only and such that the transfer function from these inputs to all other variables 
is proper (see [20^ Remark 5.11]). In view of these facts and the inequality mCB) ^ (t_|-(S), we assume 
without loss of generality 
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where m is the number of input 
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and define J, 



pq 







(4) 



3.1 Dissipativity of uncontrollable behaviors: main results 

The following theorem assumes an unmixing condition - no pair of uncontrollable poles are symmetric 
with respect to the imaginary axis. If this unmixing condition is satisfied for the uncontrollable poles, 
then the controllable part of a behavior being dissipative is equivalent to the dissipativity of the whole 
behavior. 

Theorem 3.2 Consider a behavior 53 G and a nonsingular, symmetric S G W"^'" with the input 
cardinality of ^ at-most the positive signature ofT,, i.e., m(5S) ^ cJ-|-(S). Assume that the uncontrollable 
poles are such that Aun H — Au„ = 0. Let 'B have an observable image representation, *B = Image M{-^), 
where M(^) is partitioned as 

Wi G M"""" [C] , W2 G M^P^s)^" [^] . (5) 

Let G{S,) := W^2(0^i(0""^ c-'^d D := lim^^oo G{juj). Assume ^cont is such that {Im + D'^JpqD) > 0. 
Then, 55 is T,-dissipative if and only if its controllable part ^cont is T,- dissipative. 



M{0 
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The above result gives conditions under which the autonomous part of a behavior plays no hindrance 
to dissipativity of the behavior after the controllable part is dissipative. One of the conditions for this 
is that the uncontrollable poles are not 'mixed', meaning no two of the uncontrollable poles add to zero. 
For autonomous LTI systems, this condition is a necessary and sufficient condition for solvability of the 
Lyapunov equation. Of course, storage functions are just generalizations of Lyapunov functions to non- 
autonomous systems. The other condition: (/^ + JpqD) > on the controllable part is a kind of 
strictness of dissipativity 'at the infinitj|^ frequency'. While this condition allows the use of Hamiltonian 
matrices in the proofs, this condition also rules out consideration of lossless systems from the above result. 
More significance of the assumption is noted in the remark below. 

Remark: While a behavior *B € £" admits many i/o partitions, and for each i/o partition, admits many i/s/o 
representations, the matrix / + Jq pD depends only on the behavior, and in particular, only on 'Bcont the 
controllable part. In other words, this matrix can be found from any image representation of Scont- This matrix 
being positive semi-definite is a necessary condition for dissipativity of *Bcont , and this denotes dissipativity at very 
high frequencies, i.e. as a; — >■ oo. Positive definiteness of / + D^Jq pD may hence be termed as 'strict dissipativity 
at the infinity frequency'. This assumption helps in the existence of a Hamiltonian matrix, our proofs use the 
Hamiltonian matrix properties and techniques from indefinite linear algebra (see |10|). Positive definiteness of 
this matrix is guaranteed, for example, by strict dissipativity of a behavior; on the other hand, lossless controllable 
behaviors have this matrix as zero. Another interpretation of the condition I+D^ Jq^pD > is that the 'memoryless' 
part of the behavior is strictly dissipative; the memoryless/static part of a behavior was defined in [21j, and we 
don't require this notion in this paper. 

We now review some existing results and formulate/prove new results about Hamiltonian matrices in 
the context of dissipativity of controllable and uncontrollable systems. These are required for the proof 
of the main result. 

To prove that the behavior !B is S-dissipative, we use the proposition below and show that the existence 
of a symmetric solution to the dissipation LMI is sufficient for S-dissipativity. Let the behavior have an 
input/output partition w = {wi,W2), wi G (i°°(]R,M'") as the input and W2 G (M, RP'^'i) as the output. 
Let the input/state/output representation of the behavior *B be, 

—x = Ax + Bwi, W2 = Cx + Dwi. (6) 
with {C,A) observable. The following result is well-known: see \19 \ \n \ [5]. 

Proposition 3.3 Let a behavior 5S G have an i/s/o representation for an input/output partition 
w = {wi,W2) with A,B,C,D as state space matrices and {A,B) possibly uncontrollable. The behavior *B 
is T,- dissipative if there exists a symmetric K G R"^" which solves the following LMI 

' {KA + A^K - C^JpgC) (KB - C^JpgD) 
{KB - C^JpgDf -(4 + D^JpgD) 

Further, under the condition that {A,B) is controllable, existence of K solving the above LMI is necessary 
and sufficient for T,-disspativity of^. 

^The feed-through term D is finite, since the transfer function is proper; see text before Equation Q. 



^ 0. 



(7) 
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As (Im + D^JpqD) is invertible, the Schur complement of (I^ + JpgD) in the above LMI gives the 



Algebraic Riccati Inequality (ARI) 



K{A - B{I^ + D^JpgD)-^D^JpgC) 



+{A - B{U + D^JpqDr'D^JpgCfK 



(8) 



+KB{I^ + JpqD)-^B' K-C {Jpq + DD' y^CA 
The corresponding equation is the Algebraic Riccati Equation (ARE) and we use properties of this ARE 



in proving Theorem 3.2 We rewrite the above ARE as 



where 



KA + A^K + KDK - C = 



A:={A- B{h + JpqDY'D' Jp,C) 



(9) 



and 



D := B(4 + D'^JpgD)-^B^, C := C'^ {Jpq + DD^Y^C. 
Also define the Hamiltonian matrix corresponding to the ARE: 



H 



A D 

C -A* 



(10) 



We define the following matrices; they play a crucial role in the results we use from ^lOj and in our 
proofs. 



M := j H, P :-- 



and P := j 



/ 
-/ 



(11) 



-C A* 
A D 

It is well-known that a symmetric solution to ARE can be obtained from an n- dimensional, M- 
invariant, P- neutral subspace and we state this in the proposition below. Let K G be a Hermitian 

matrix. The graph subspace corresponding to matrix K is defined as, 



SiK) := Im 



I 
K 



Notice that '5{K) is an n-dimensional subspace of C^". The following proposition from [lOj states solv- 
ability of the ARE in terms of P-neutrality with respect to the above P. 

Proposition 3.4 (See TWp Consider the ARE ^ with {A,b) possibly uncontrollable. Then K € C"^" 
is a Hermitian solution of the ARE if and only if the graph subspace of K is M -invariant and P -neutral. 



The proposition below allows use of simplified A and B for all later purposes when dealing with 
uncontrollable systems. 

Proposition 3.5 (See ]2^) Consider the behavior OS G ii"' with an input/ state/ output representation 
= Ax + Bwi, W2 = Cx + Dwi, where w = {wi, W2). Then there exists a nonsingular matrix T S M"^" 
such that 

Cc Cu 



T-^AT 



Ac 




Acp 
A-u 



T"'B 



Be 




and CT 
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Further, 

— X = AcX + BcWl, W2 = CcX + Dwi. 
at 

gives an i/s/o representation for the controllable part ^cont- 



Let He and Mc denote the corresponding matrices for the controllable part ^cont as defined in (11) 



with Ac,Dc and Cc defined accordingly in (10). 



Lemma 3.6 Suppose ^cont G ^cont satisfies the assumption that {Im + D^JpgD) > 0. // ^cont is S- 
dissipative, then the partial multiplicities corresponding to the real eigenvalues of Ale, if o-ny, are all even. 



In order to prove Lemma 3.6, we use a result from [lOJ concerning the partial multiplicities of real 
eigenvalues of M. Let the generalized eigenspace of a matrix A corresponding to an eigenvalue Aq be 
denoted by Jlx^^,{A) and the controllable subspace of the pair {A, D) be denoted by ^. 

Proposition 3.7 (See JJEI) Consider the behavior ^cont S ^cont- Assume Dc > and C* = Cc and 

there exists a Hermitian solution K E C^^"^ to the ARE Suppose 

3ix,{Ac + DcK)0^^^j^^ 

for every purely imaginary eigenvalue Aq of {Ac + DcK). Then, the partial multiplicities of corresponding 
real eigenvalues of Mc are all even and are twice the partial multiplicities of those corresponding to the 
purely imaginary eigenvalues of (Ac + DcK) . 



Proof of the Lemma 3.6 As the controllable part is S-dissipative, there exists a symmetric solution 



K to the ARE by Proposition 3.3 Since the behavior is controllable, the following is true: ^ = C" 
and hence Jlx^{Ac + DcK) C for every purely imaginary eigenvalue Aq of {Ac + DcK). Thus, using 

Proposition |3.7[ the partial multiplicities of all real eigenvalues of Mc, if any, are all even. This completes 
the proof. 

We define a set called c-set as in [10]. Such a c-set, if exists, guarantees the existence of a unique 
P-neutral, M-invariant subspace N under certain conditions. This is made precise in the proposition 
below. 



Definition 3.8 (See 1 10^ ) Let M £ {[^nxn Q be a finite set of non-real complex numbers. Q is 

called a c-set of M if it satisfies the following properties 

1. ene = 

2. C U C = cr(M) \ M, the set of all non-real eigenvalues of M . 

Let be an invariant subspace of M and denote the restriction oi A to Tsf as A\yi. 

Proposition 3.9 (See JTUj) Let M G ^'irix2n ^ P -self-adjoint matrix such that the sizes of the Jordan 
blocks of M , say: mi,m2, ■..,rnr, corresponding to real eigenvalues of M are all even. Then for every 
c-set C there exists a unique P-neutral M-invariant subspace of dimension n and a{M\yi) \]R = C, and 
the sizes of the Jordan blocks of M\yi corresponding to the real eigenvalues are ^ mi, ^ m2, ^ m^. 
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It is easy to verify that M and P satisfy the fohowing relation, PM = M*P. Hence M is P-self- adjoint 
and Proposition |3.9| can be used. 

Following proposition, which is a reformulation and combination of Theorems A. 6.1, A. 6. 2 and A. 6. 3 
in [10], states that the partial multiplicities of an eigenvalue are unaffected by pre- multiplying and/or 
post-multiplying by an unimodular matrix. This result is used in the proof of the Theorem |3.2[ 



Proposition 3.10 Consider Si{(,) and S'2(0 e and let and pj{0 £ for i = 1, . . . ,w 

be the invariant polynomials of Si{^) and S2{6.) respectively. Suppose Si{^) = TiS2{£,)T2, for invertible 
Ti, T2 G M"^™. Let A G C and . . . , be the maximum integers, for j = 1,2, such that — X)^^ divides 
P)iC) fof j = ^, ■ ■ ■ ,u!- Then I3j = (3j for j = 1, . . . ,w. In particular, if 



and 



ii-p 

Q{i) 

^I-P 
I 



and det(5(A) 7^ 0, then the partial multiplicities of Si{^) and 52(^) corresponding to X are equal. 



Before proving the main result Theorem 3.2, we state and prove another useful result. The follow- 
ing lemma relates the partial multiplicities of purely imaginary eigenvalues of the Hamiltonian matrix 
corresponding to the controllable part to that of the uncontrollable behavior. 



Lemma 3.11 Consider the behavior 55 G with the set of uncontrollable modes A„„ satisfying A„„ n 
—Aun = 0- Let 55 have an observable i/s/o representation as = Ax + Bwi, W2 = Cx + Dwi, 
induced by w = {wi,W2), such that {Im + D^JpqD) > 0. Further, let ^cont = ^™^(^)- Define 
^{CiV) — M'^ {(^)TiM{ri) and construct the Hamiltonian matrix, H as in (10). Then 

1. a{H) = roots (det U A„„ U (-A„„) . 

2. If the controllable part ^cont is T,-dissipative, then the partial multiplicities corresponding to the 
purely imaginary eigenvalues of H , if any, are all even. 



Proof of Lemma 3.11 Statement 1: See 
Statement 2: We use the fact that if the controllable part is dissipative, then the partial multiplicities 
corresponding to purely imaginary eigenvalues of He are even. Without loss of generality the following 
i/s/o representation for 5S is assumed 



d^ 
dt 









A 




Xq 


+ 





































W2 



Cc Cu 



with [Ac, Be) controllable. 

Then, the Hamiltonian matrix gets the following form 



Xc 
Xu 



+ Dwi 



Wl 



(12) 



(13) 
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H 



Ar 



A. 



cp 











Cj^Cc c'^Cu ~a'^ 

fiTfi (-,Tfi _aT _aT 



where 



Acp — A(yp ^c(-/ni ~(~ Jpq-b^ b Jpq^ixi 



and 



Be = B,{I^ + D^JpgD) 



Cc = (/pq + DD^yWe, C„ = (Ipq + DD'Y^^Cu, 



and Au = Au- 

The Hamiltonian matrix for Q3cont is 



CjCc -AT 



Consider the polynomial matrix H{^) = ^l2n — H, 

^Ific ~ Ac —A, 



cp ~BcB^ 
^Inu-Au 

—cfTCr —cJ!Ci,, aJ ^Iuu + 



Applying the following transformations to H{^), we get -ffi(C) = EiH{^)Ei 



ic Ac B^B^ Acp 
CcCc Clfic + Ac —CjCu 

- v4„ 









AT 
cp 



~C^Cu ^luu + A^ 



where 



El :-- 



In,, 

In, 

/n„ 

In., 



(14) 



Now, since A„„njM = 0, for A G a{Hc)r\j'M., the matrix blocks A/„^ — Au and XIuu ^^^^ invertible. 
Pre-multiplying -ffi(0 by -E2 and post-multiplying by £^3, we get -^2(0 '■= B2Hi(^)E^ = 



iln, - Ac -BcBj 

Cln^ - Au 


















<Cu iln^ + Al 



(15) 
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where 



In. 


"J 1 





In. 










Inu 











In 



and 



In. 














In. 














Inu 












Inu 



and Ti := A/^^ - A„, r2 := A/„„ + , Ts := -Acp and T4, := -C^^ C„. 
Thus 

' e^2n, - 



H2{0 = E2EiH{0EiE2 







QUO 



where 



Qu 



Ti 



Thus by using Proposition 3.10, the partial multiplicities of purely imaginary eigenvalues of H are 



even. This completes the proof of Lemma 3.11 



4 Proof of the main result: Theorem 3.2 



Proof: If part: Assume that the controllable part ^cont is S-dissipative and the assumption in the 



theorem is satisfied. By Propositions 3.3 and 3.4, to prove that the behavior !B is dissipative, it suffices 
to show the existence of a G C"^° such that the corresponding graph subspace is an n-dimensional, 



M- invariant, P- neutral subspace. To show the existence of such a K, we use Lemma 3.9 to construct a 
c-set such that the corresponding n-dimensional M- invariant, P- neutral subspace is also a graph subspace 
of K. This K is the solution to the ARE and storage function for the whole behavior would then be 
defined Kx, thus completing the proof. 

As the unmixing assumption on uncontrollable modes is assumed, A G jMr]a{H) means that A ^ A^^ 



and A G a{Hc)- As we have assumed that the controllable part is S-dissipative, from Lemma 3.11 the 



partial multiplicities of real eigenvalues of M{:= iH) are all even. Using this fact, Lemma 3.9 can be used 



to infer that there exists a unique n-dimensional M- invariant, P- neutral subspace for every c-set. 

Now, it remains to show the existence of a c-set such that the corresponding n-dimensional, M- 
invariant, P-neutral subspace is also a graph subspace. 

We choose a c-set C such that jAu„ C C and show that the corresponding n-dimensional, M-invariant, 
P-neutral subspace is a graph subspace. Let L be the n-dimensional, P-neutral, M-invariant subspace of 
C^" corresponding to the c-set S and suppose 



L = Im 



X2 



(16) 
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for matrices Xi and X2 G C°^°. In order to prove that £ is a graph subspace it is enough to prove that 
Xi is invertible. This is proved using contradiction: we assume Xi is singular and show that we get a 
contradiction to the unmixing assumption on Kun- This constitutes the rest of the proof of the 'if part'. 
We prove this along the lines of |10j . 
M-invariance of L implies that 



" A 


D 




' Xi ' 




" Xi ' 


C 


-A* 




. ^2 




. ^2 . 



for T G C^^". In other words, 



Then as L is P-neutral 



j{AXi + DX2) = XiT, j{CXi - A*X2) = X2T. 



X^ X2* 



' -C A* ' 




' Xi ' 


A D 




. ^2 . 



(17) 



(18) 



i.e. X^DX2 + XIA*X2 + X^AXi - X^CXi = 0. 



Now suppose Xi is singular. Let % = kerXi. From (18), for every x E % we have 



x*X^DX2X + x*XlA*X2X + x*X^AXix - x*XlCXix = 0, 



which implies 



x*XWX2X = 0. 



(19) 



Since D >0, X2X G kerD, i.e., X2% C kerD. Now, for every x ^ %, from equations (17) we have, 

XiTx = jAXix + jDX2X = 0, 

that is, 

TX C %. 



(20) 



This implies % is T-invariant. Hence there exists a non-zero v which is an eigenvector of T such that 
Xiv = corresponding to eigenvalue, say A. We claim that A cannot be a real eigenvalue and prove this 
below. 



Post-multiplying the second equation of (17) by f we get 



jCXiv - jA*X2V = X2TV. 
- jA*X2V = XX2V. 



(21) 
(22) 



This implies that X2V is a left eigenvector of A with eigenvalue —j A and from equation (19) we have 
B'^X2V = 0. This means that —j A is an uncontrollable eigenvalue of A and —j A G i-G-j A G jA^^. 
Now, if A were a real eigenvalue, then —jX and j A belong to Aun and contradicts the unmixing assumption 

Kn n Aun = 0. 

Since o"(T) \M = C (by Proposition 3.9), A G C. We now have A G C but A cannot belong to C. But we 
have A G jAun and thus A G C which is a contradiction by definition of c-set. This completes the proof 
for the if part. 
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Only if part: Assume ^ is S-dissipative. Then there exists a storage function Q^{w) such that 

d 



dt 



Q^{w) ^ Qt.{w) 



Integrating both sides for every ti; G 55 n !D , 



QT,{w)dt ^ 



which implies that 5Scont is S-dissipative. This completes the proof of Theorem 3.2 



(23) 

(24) 
□ 



The above proof is constructive in the sense that if a behavior !B G satisfies the three conditions: 

• uncontrollable poles are unmixed, i.e. no two of them add to zero 

• the controllable part !Bcont is dissipative, 

• the controllable part QScont is 'strictly dissipative' at infinity, i.e. (1^ + D^JpqD) > where D is 
the feed-through term of the transfer function, 

then we construct a storage function that satisfies the dissipation inequality for the whole behavior !B. 
Further, the storage function we construct is equal to Kx where X is a solution to the correspond- 
ing Algebraic Riccati Equation. Further, we constructed the storage function by starting with a state 
representation in which the state x is observable from the manifest variable w. These facts lead to the 
following important corollary. 

Corollary 4.1 Let *B G 6e an uncontrollable behavior whose uncontrollable poles A„„ are unmixed, 
i.e. Kun n — Au„ = 0. Consider S partitioned in accordance with the input cardinality of ^ as 

'L 



Ig 









Suppose *Bcont, the controllable part of^, has an observable image representation, w = M{-^)£, where 
M(^) is partitioned as 

WiiO 



M{0 



Wi G 



'[i\,w2 G 



o(y+g)xm 



Let G{s) := W2{s)Wi{s) ^ and D := lim^-^oo G(s). Assume ^cont is such that {Lm + JpqD) > 0. 
Then, the following are equivalent. 

1. 53cont is dissipative. 

2. There exists a ^' G M™^"'[C,'?/] such that Q^{w) is a storage function, i.e. ^Q-^{w) ^ w'^TiW for all 
t/; G 53. 

3. There exists a matrix K and an observable state variable x such that ^x^Kx ^ w'^T.w for all 

Statement 2 tells that the storage function can be expressed as a quadratic function of the manifest 
variables w and their derivatives. Statement 3 says that the storage function is a 'state function', i.e. a 
static function of the states, and hence storage of energy requires no more memory of past evolution of 
trajectories than required for arbitrary concatenation of any two system trajectories. 
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5 Examples 



In this section we discuss two examples of uncontrollable systems that are dissipative. The Riccati 
equations encountered in these cases are solvable by the methods proposed in this paper; we also give 
solutions to the Riccati equations. 

The first example is of an uncontrollable system with uncontrollable modes satisfying the unmixing 
assumption, i.e. no two of the uncontrollable poles add to zero. However, the Hamiltonian matrix has 
eigenvalues on the imaginary axis. 

Example 5.1 Consider the behavior *B whose input/ state/ output representation is given by the following 
A, B, C and D matrices 



A 






-0.5 




-0.5 






, B = 




1 


-1.5 




-0.5 



, c 



-0.5 



D = 0.5 



with (t{A) = { — |, —1}- Here K^n = {—1} which satisfies the unmixing assumption. An equivalent kernel 
representation of the behavior is given by 







(^2 + 2^ + 1) -(2^2 + 3^ + 1) 
In this case S = diag(l,— 1) and hence = m(?B), it can be checked that the controllable part 



Wl 
W2 



cont 



ker 



(e + 1) -(2^ + 1] 



3.2 



is Tj-dissipative. And [Im + D) = 3/4 > 0. Thus from Theorem 

!B is Tj-dissipative. 

The following real symmetric matrix induces a storage function that satisfies the dissipation inequality 



K 



1 



7 -1 
-1 1 



im 



The 2- dimensional, M -invariant, P-neutral subspace which gives the solution is 

j 8 
3 1 
1 

The next example is an RLC circuit shown in Figure [TJ 

Example 5.2 Consider the RLC circuit system whose input u is the current flowing into the circuit, and 
output y is the current through the inductor. A state space representation of the system is found using 
the following definition of the states. The state variables are xc, voltage across the capacitor, andxi, the 



ductor. Assume Rl 


= Rc =: 


R. 




A = 


' -1/c" 


, B = 


'l/C 


, c 




1/L -2R/L 




R/L 





1 



D = R, 



The system becomes uncontrollable when L = R^C. Let R = 0.5, C = 1 and L = 0.25. The state 
representation of the system is 
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Figure 1: An RLC circuit system with input u and output y 



A 



-1 




1 








4 -4 




2 



c 



1 



i:' = 0, 



Poles of the system are —2, —2 and one o/ them is uncontrollable. The controllable part is dissipative 
and the corresponding Hamiltonian matrix has eigenvalues on the imaginary axis. A solution to the ARE 
is given by the following symmetric matrix 



K 



3 -0.5 
-0.5 0.25 



(25) 



which induces the storage function — xcxl + x\/A. 



6 Behaviors with static controllable part 

In this section we consider the case when all the states of the system are uncontrollable, in other words, 
when the only controllable part is static/memory less. In a state space realization, this means that the 
matrix B is zero. The case of autonomous behaviors is clearly a special case: D also is zero and the 
assumptions in the lemma are satisfied. 

Lemma 6.1 Consider a behavior *8 with static controllable part. Let the behavior have a state represen- 
tation, 

—x = Ax, W2 = Cx + Dwi. (26) 

with the pair (C, A) observable. Assume {Im + D^JpqD) > and A„„ C iM.. Then there does not exist a 
symmetric solution to the corresponding ARE. 

Proof: The Hamiltonian matrix takes the form, 

C -A* 

Here A = A, C = C^{Jpg + DD'^Y^C. 
Next we use the following proposition to say that partial multiplicities of purely imaginary eigenvalues 
of the Hamiltonian matrix are even. 



(27) 
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Proposition 6.2 Consider the matrix 



N 



A 




D 

-A* 



(28) 



Then for every purely imaginary Aq G c(^) such that 3lx^^{N) C ^^f), the partial multiplicities of such 
Ao are even. In fact, they are twice the partial multiplicities of Aq as an eigenvalue of A. 



Consider, 



A^ 

= . . (29 

-a\ 

As the pair {C,A) is observable, the pair {A^,C^) is controllable. From Proposition 6.2, the partial 
multiplicities of purely imaginary eigenvalues are twice the partial multiplicities of purely imaginary 
eigenvalues of A. 



7 Lossless autonomous behaviors 

In this section we investigate the requirement of unobservable variables in the definition of the storage 
function. As has been studied/shown so far, for controllable dissipative systems (|20j). the storage function 
need not depend on unobservable variables. It was later shown in [12] that for the case of strict dissipativity 
of uncontrollable systems, observable storage functions are enough under unmixing and maximum input 
cardinality conditions. Theorem |3.2| shows that this is true for a more general scenario, i.e., dissipativity 
(including non-strict dissipativity) for all input cardinality conditions under unmixing assumption. On the 
other hand, when relaxing the unmixing assumption elsewhere except on the imaginary axis, under certain 
conditions solutions to the ARI exists though ARE does not have a solution (see [9]). In this section, we 
investigate the need for unobservable storage functions for uncontrollable systems whose uncontrollable 
poles lie entirely on the imaginary axis. We discuss the case for autonomous behaviors below. 

Lemma 7.1 Consider an autonomous behavior ^aut 

—X = Ax, w = Cx, 
at 

with cr{A) Pi iM 7^ 0. Let the supply rate he Qt.{w) = —w^w. Then the following is true. 
If there exists a storage function Q^{w) satisfying the inequality 

j^Q^iw) i^Qj^iw), ywe'Baut (30) 
then any A G cr(yl) n zM is C-unobservable. 

Proof: Suppose if there exists a storage function which is a state function x'^Kx satisfying the dissipation 
LMI, then the dissipation LMI ([T]) is equivalent to the Lyapunov inequality 

KA + A*K + C*C i^O. (31) 

Now for every eigenvector of A corresponding to eigenvalue A G zM, we have 

x*iKA + A*K + C*C)x s^O, (32) 
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which gives 

Xx*Kx + Xx*Kx + x*C*Cx(^0 (33) 

or 

x*C*Cx ^ 0. (34) 

But, as x*C*Cx > 0, we have Cx = for every eigenvector of A corresponding to A G cr{A) D iM. This 
implies that the any A G cr{A) (1 iM is C-unobservable. 

This observation teUs that for dissipativity of autonomous systems having eigenvalues on the imaginary 
axis, it is necessary to allow storage functions to depend on unobservable variables also. 

8 Orthogonality and uncontrollable behaviors 

In this section we investigate the property of orthogonality of two behaviors in the absence of controlla- 
bility. We propose a definition that is intuitively expected and show that by relating this definition to 
lossless uncontrollable behaviors, we encounter a situation that suggests an exploration whether dissipa- 
tivity should be defined for behaviors for which the input-cardinality condition is not satisfied. We first 
review a result about orthogonality of controllable behaviors. 

Proposition 8.1 Let S G ]^«'X"' nonsingular, and suppose *Bi,*B2 G -Ccont- The following are equiv- 
alent. 



1. wlT.W2dt = for all wi G^iDD and for all t(;2 G *B2 n D. 

2. 5Si X *B2 is lossless with respect to 



S 




3. There exists a bilinear differential form L^, induced by ^ G M'"^"'[C,??] such that ^L^{wi,W2) = 
Statement 1 above is taken as the definition of orthogonality between two controllable behaviors !Bi 



and *B2 in [20]. Keeping in line with Definition 3.1 for dissipativity, we could take Statement 3 above as 
the definition of orthogonality for behaviors not necessarily controllable. The drawback of this approach 
is elaborated later below in this section. We pursue a different direction as follows. Notice that if !Bi 
and 532 satisfy the integral condition in Statement 1, then this integral condition is satisfied for every 
respective sub-behaviors *B'^ and also. Of course, restricting to compactly supported trajectories in the 
integration implies only controllable parts of respectively ^B'^ and satisfy orthogonality. The following 
definition builds on this property. 

Definition 8.2 Consider a nonsingular E G M'"^'^ and let *Bi and 532 £ -2". Behaviors 53i and 532 c^re 
said to be Ti-orthogonal (and denoted by 53i _Ly; 532^^ if there exist 53f and ^B^ G -C"cont such that 

• 5Si C 53f, 

• 532 C 53^, and 

• /jj w'[T.W2dt = for all wi GRID'S and for all W2 G ^2'^'S. 
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While this definition is not existential in the storage function, it is existential in 5Sf and !B2, raising 
questions about how to check orthogonality. Note that if 55 is an uncontrollable behavior, then any 
controhable 55^= e £"cont such that 55 C QS'^ satisfies ni(5S) < ni(5S'=), and 55 ^ 

The question arises as to how much larger a controllable 55^^ would have to be for it to contain !B. 
This problem is addressed in the following subsection. 



8.1 Smallest controllable superbehavior 

Due to its significance for determining whether two uncontrollable behaviors are orthogonal, in this 
subsection we study the following problem: 

Problem 8.3 Let *Bi € Find ^2 G -2" such that 

1. *B2 5 *Bi 

2. 5S2 G -Ccont i-G., a controllable behavior 

3. is a behavior with the smallest input cardinality satisfying Properties 1 and 2. 
The following theorem answers this question. 

Theorem 8.4 For S G the following statements are true. 



1. There exists 5^2 G -C"cont satisfying the requirements in Problem 8.3. 

2. The behavior ^2 is unique if and only i/5Si is controllable, and in that case 'Bi = 5^2. 

3. Assume 5Si is uncontrollable. The input cardinality 0/ 5^2, ni('B2) satisfies m(5S2) > m(5Si). More 

precisely, m(5S2) = ni(*Bi) + k where k := max rank (i?i(A)) — min rank (i?i(A)) where Ri is a kernel 

agc Aec 

representation 0/ 5Si . 

Proof. (1): This is shown by constructing such a ^2- Let keri2i(,^) and keri?2(C) be the kernel 
representations of 5Si and 5^2 respectively where Ri G RPi^"", R2 G RpsX"" and there exists an F G RPaxpi 
such that FRi = i?2- Let Au„ be the set of uncontrollable modes of !Bi. Without loss of generality, we 
assume Ri to be 



Ri 



S 



(35) 



where S is of the form 



such that / G r(p-'^)^(p-^) is identity matrix and D G M^^''[^] is a 



/ 
D 

diagonal matrix with di,d2 . . ■ ,d\ along its diagonal and satisfying the divisibility property: di\d2, d2\d-i, 
. . . and (ik-i|c^k, with degree of di at least one, and with 



k = max rank (i?i(A)) — min rank (i?i(A)). 



Partitioning F 



Fi F2 



conforming to the row partition of Ri , we have 



R2 



Fi Fo 



10 
D 



(36) 
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This simplifies to 



i?2 



Fi F2D 



(37) 



Let be the behaviour defined by the kernel representation of R2. For 5^2 to be controllable, i?2(A) 
needs to have full row rank for all A G C. As F2P loses rank for A G A„„, Fi should have full row rank 
for every A € C so that 5S2 is controllable. This proves the existence of 5^2 satisfying properties 1 and 2 



of 8.3 In order to satisfy property 3 in Problem 8.3 i.e. m(5S2) has to be the least, we have to choose a 



unimodular Fi and free F2 such that 5S2 satisfies the three conditions. F2 can be freely chosen because 
the choice of F2 does not affect the input cardinality of 532- 

(2) : Let 5Si be controllable. Then, in the above, S = I and D does not exist. This means that F2 does 
not exist. Thus the kernel representation matrix of ^2 would be given by [-Fi 0] where Fi is unimodular. 
Therefore, 5^2 is unique. Let QSi be uncontrollable. Then D exists and as shown above, free F2 can be 
chosen. This makes the behaviour 5^2 non-unique. 

(3) : The input cardinality of 5^2 is given by 

m(*B2) = w — p2 

= w-{pi-^) (38) 
= m(QSi) +k 

Here, k = max rank (i?i(A)) — min rank (i?i(A)) which is determined by the size of the D matrix. □ 



8.2 Superbehaviors and orthogonality 

We saw above in this section that orthogonality of two uncontrollable behaviors is defined by requiring 
these uncontrollable behaviors to be sub-behaviors of two orthogonal controllable behaviors. Using the 
result on existence of super-behaviors that are controllable, and their non-uniqueness even if the are the 
smallest controllable superbehavior, we formulate the question of whether two uncontrollable behaviors 
are orthogonal as a question of finding a pair of smallest controllable superbehaviors that are mutually 
orthogonal. The requirement of them being smallest is motivated by the fact that orthogonality of two 
controllable behaviors imposes an upper bound on their input cardinalities: this is reviewed below. For a 
behavior 5S G £" and a nonsingular matrix S, the set S5S is defined as follows 

S*B := {w G e:°°(M,M^) | there exists 7; G OS such that w = Sv}. 

It is straightforward that is also a behavior, its controllability is equivalent to that of 2S, and the 
input cardinalities are equal. 

Proposition 8.5 Let QSi and 052 £ -2"cont CLi^-d suppose E G R"'^'^ is nonsingular. Then, the following 
are true. 

1. *Bi _Le *B2 4^ «Bi _L (SQS2). 

2. *Bi _L 5S2 ^ m(<Bi) -Mn(5S2) ^ w. 
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3. *Bi ±E *B2 ^ ni(QSi) +m(«B2) ^ w. 



Due to the above inequality constraint on the input-cardinahties of orthogonal controllable behaviors, 
the uncontrollable behaviors too have a necessary condition to satisfy for mutual orthogonality. 

Lemma 8.6 Suppose *Bi and 5^2 G with at least one of them uncontrollable and let S G M"'^"' he 
nonsingular. Assume *Bi J-s *B2- Then ni(5Si) +ni(5S2) < w. 

Example 8.7 Consider the pair of 'seemingly' orthogonal behaviors studied in 12^. page 360]. Define 5Si 
and ^2 G by 

= Ax, wi = Cx, W2 free, i.e. B2 = (J:~(M,M") 

and the supply rate wj'w2. Thus 5Si is autonomous. Consider the following non-observable latent variable 
representation for ^2-' gj-z = —A^z + C'^W2. It can be checked that the 'storage function' x'^z satisfies 
z = W1W2. Existence of such a storage function is, in fact, reasonable for a (different) definition 
of orthogonality of the two behaviors QSi and 552- In fact, any autonomous behavior *Bi G £" is then 
'orthogonal' to 5^2 = C°°(M,M^)/ However, the 'embeddability' definition we have used above rules out 
this example for an orthogonal pair of behaviors since the necessary condition of the above lemma is not 
satisfied. In other words, there doesn't exist a controllable behavior such that contains *Bi and 
^1 J-S ^2- Thus 5Si and ^2 o-fe not I -orthogonal. □ 



9 Dissipative sub-behaviors/superbehaviors 

In this section we look into the input cardinality condition for dissipative behaviors. Recall that a behavior 
*B G which is dissipative with respect to the supply rate S (constant, symmetric, nonsingular matrix) 
satisfies the condition iii(*B) ^ cr+(S). We now look into the possibility of embedding a behavior OS in a 
controllable superbehavior that is S-dissipative, and into the drawback of using this as the definition of 
dissipativity, along the lines of orthogonality defined in the previous section. 

Problem 9.1 Civen a nonsingular, symmetric and indefinite S G M'^^"', find conditions for existence of 
a behavior 'B G £" such that 

• there exist and !B_ G -C"cont 

with *B = «B+ n *B_. 

• 53+ is strictly S dissipative 

• 53- is strictly — S dissipative. 

The significance of the above problem is that if a nonzero behavior 53 satisfying above conditions exists, 
then clearly such a behavior would be both strictly S and strictly — S dissipative, raising concerns about 
whether embeddability in a dissipative controllable superbehavior is a reasonable definition of dissipativity 
(when dealing with uncontrollable behaviors). The following theorem states that nonzero autonomous 
behaviors can indeed exist satisfying above condition. 

Theorem 9.2 Let S G M'"^"' be nonsingular, symmetric and indefinite. Then there exists a nonzero 



!B G -C" such that requirements in Problem 9.1 are satisfied. Further, any such !B satisfies m(!B) = 0, i.e. 
!B is autonomous. 
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Proof. We first show that given any nonsingular, symmetric and indefinite S G M"'^"', a behavior OS 
satisfying above properties exists. Without loss of generahty, let 









-/_ 



with sizes of the identity matrices /+ and /_ equal to and a^. Due to indefiniteness of S, 0"+, (T_ ^ 1 
and, due to the nonsingularity, they sum up to w. Choose anjj^ nonzero polynomial matrix G 
M^^™+[^] with Ml{-joj)T.M+{iuj) ^ e+/^ for all G M for some e+ > 0. Similarly, choose M_ G 
M"'^™-[^] such that MT(-ja;)i;M„(ja;) ^ -e_/^„ for all a; G M for some e_ > 0. Now define 03+ and 
!B_ by image representations w = M+(gj)£ and w = M^{^)£ respectively. Define 5S := 53+ n 53_, thus 
proving existence of 53 as stated in the theorem. 

We now show that !B is autonomous. Let i?+ and R- be minimal kernel representation matrices of 
and !B_ respectively. Then !B = !B+ n 5S_ is described by the kernel representation matrix R G M*^"'[.^] 
R 

Clearly, rank (R) ^ w. Suppose rank (R) < w. Then, there exists p G 



with R := 



R. 



and 



p ^ such that i?+p = and R-p = 0. This implies Image G and Image G 5S_. 

Taking w = with £ G D(M, M) and £ / 0, it follows that w G D(M^, M). Further, p / 0, hence t/; / 

because ^ is non-zero and of compact support. Further, we have e+, e_ > such that 



w 



+ \M\L2 



and 



w 



kill 



Both the above conditions cannot be satisfied simultaneously for w ^ Q. Thus, rank [R) < w gives a 
contradiction. This proves rank (R) = w and hence autonomy of !B. □ 



We illustrate the above theorem using an example. 

Example 9.3 Let S = diag (1,-1). Define 53+ and !B_ by image representations w = M+(gj)^ and 
w = M-{47)£ respectively, with 



e + 4 

3 



and M_ (^) 



2 

e + 5 



Strict dissipativities is easily verified. Calculating the kernel representations, we get a kernel representation 
for 53 := «B+ n 53_ as R{^)w = 

-2^ 

Clearly, R is nonsingular and hence 53 is autonomous. □ 



For non-strict dissipativity case, we have the following problem and theorem. 



^Such matrices are plenty due to existence of sufficiently many controfiable strictly dissipative behaviors for every supply 
rate E satisfying the stated conditions. 

^It is not difficult to show that if M+ and M_ were nonconstant polynomial matrices, then 2J is not the zero behavior. 
An example following the proof makes this easier to see. 
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Problem 9.4 Given a nonsingular, symmetric and indefinite S G M"'^"', find conditions for existence of 
a behavior !B G £" such that 

• there exist 53+ and !B- G Xl^cont suc/i that <B = <B+ n <B_. 

• 53+ is S dissipative 

• 53- is -S dissipative. 

Theorem 9.5 Zei S G M"'^"' 6e nonsingular, symmetric and indefinite. Then there exists 53 G £" such 



that requirements in Problem 9.4 are satisfied. Any such 53 satisfies m(5S) ^ min(cr+(S), cr_(S)). In case 

53 is uncontrollable, ni(53) < min(cr+(S), cj_(S)) 

If m{^) ^ 1, i/ien neither nor^^ can be strictly dissipative. 

Proof. The proof proceeds in the same way as the proof for the previous theorem, except for the 
strictness of the dissipativities. Construct 53+ and *B_ as in the previous proof, but with e+ and e_ equal 
to zero. We have 

/oo 
uFT^wdt ^ for ah u; G 53_ n 2) 

The above two equations imply w'^Tiwdt = for all w G 'B n 5). Since 53 = ^3+ n *B_, the behavior 
*B is dissipative with respect to both S and — S. Dissipativity with respect to S implies m(!B) ^ (J+($]). 
Similarly, dissipativity with respect to — S implies m(!B) ^ o"_(S). This implies 

m(53) ^ min(CT+,cr_) (39) 



If 53 is uncontrollable, then from Theorem 8.4, then the two inequalities leading to the inequality (39) 
are both strict. Hence, the input cardinality of 53 has to be strictly less than that of 53+ and *B_. This 
implies 

m(53) < min((T+,cr_) (40) 



If m(5S) ^ 1, then from Theorem 9.2 *B+ and !B_ cannot be strictly dissipative with respect to S 
and — S respectively. This completes the proof. □ 



As one of the consequences of the above theorem, if the input cardinality condition is satisfied for an 
uncontrollable behavior, i.e. m{^) = (T+(S) or m(53) = (t_(S), then such a behavior cannot be embedded 
into both a S-dissipative controllable behavior and a — S-dissipative controllable behavior. However, an 
observable storage function for such a situation exists when the controllable part is strictly dissipative 



at oo and when the uncontrollable modes satisfy the unmixing condition (see 3.2 and |16| [T2] for this 
situation in presence of more assumptions), A situation when m(5S) = o"+(S) is very familiar: we deal 
with RLC circuits in the next section. 

6. We use this method of defining orthogonality of two behaviors to explore further the definition of 
dissipativity of a behavior. Here we bring out a fundamental significance of the so-called input-cardinality 
condition: the condition that the number of inputs to the system is equal to the positive signature of the 
matrix that induces the power supply. We show that when this condition is not satisfied, then a behavior 
could be both supplying and absorbing net power, and is still not lossless. 
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10 RLC realizability 



In this brief section we revisit a classical result: a rational transfer function matrix being positive real is a 
necessary and sufficient condition for that transfer matrix to be realizable using only resistors, capacitors 
and inductors (see [6] and also [1] for the case with transformers) . Note that the transfer matrix captures 
only the controllable part of the behavior and positive realness of the transfer matrix is nothing but 
dissipativity with respect to the supply rate v^i. This is made precise below. 

Consider an n-port electrical network (with each port having two terminals) and the variable w = {v, i), 
where v is the vector of voltages across the n-ports and i is the vector of currents through these ports, 
with the convention that v'^i is the power flowing into the network. Behaviors that are dissipative with 
respect to the supply rate v'^i are also callecQ 'passive'. Given a controllable behavior 5S whose transfer 
matrix G with respect to a speciflc input/output partition, say current i is the input and voltage v is 
the output, is positive real, one can check that this behavior is passive. G in this case is the impedance 
matrix and is square, i.e. the number of inputs is equal to the number of outputs. One can introduce 
additional laws that the variables need to satisfy, thus resulting in a sub-behavior ^suh ^ ^ which has 
a lesser number of inputs; consider for example these additional laws as putting certain currents equal 
to zero: due to opening of certain ports. However, the transfer matrix for 5Ssub with respect to the 
input/output partition: input as currents through the non-open ports and output as the voltages across 
all the ports, is clearly not square, and in fact, tall, i.e. has strictly more rows than columns. Let Gsub 
denote this transfer function. Since ^sub C 55, the behavior ^suh is also passive. Of course, 5Ssub need not 
be controllable, even if 'B is assumed to be controllable. As an extreme case, suppose all the currents are 
equal to zero, and we obtain an autonomous ^snh- While RLC realization of such transfer matrices which 
are tall, and further of autonomous behaviors obtained by, for example, opening all ports, has received 
hardly any attention, we remark here one very well-studied sub-behavior of every passive behavior: the 
zero behavior. 

Consider the single port network with v = and i = 0. This port which is called a nuUator (also 
nullor, in some literature) behaves as both the open circuit and shorted circuit: see [21 page 75] and [8]. 
The significant fact about a nullator is that a nullator cannot be realized using only passive elements, and 
moreover any realization (necessarily active) leads to the realization of both a nullator and its companion, 
the norator; a norator is a two-terminal port that allows both the voltage across it and current through 
it to be arbitrary. 

11 Concluding remarks 

We briefiy review the main results in this paper. We first used the existence of an observable storage 
function as the definition of a system's dissipativity and proved that the dissipativities of a behavior and 
its controllable part are equivalent assuming the uncontrollable poles are unmixed and the dissipativity 
at infinity frequency is strict. This result's proof involved new results in the solvability of ARE and used 
indefinite linear algebra results. 

We showed that for lossless autonomous dissipative systems, the storage function cannot be observable, 

*It is also common to require the storage function in the context of this dissipativity to satisfy non-negativity. The 
sign-definiteness of the storage function is not the focus of this paper, hence we ignore this aspect. 
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thus motivating the need for unobservable storage functions. We then studied orthogonahty/lossless 
behaviors in the context of using the definition of existence of a controllable dissipative superbehavior 
as a definition of dissipativity. In addition to results about the smallest controllable superbehavior, we 
showed necessary conditions on the number of inputs for embeddability of lossless/orthogonal behaviors 
in larger controllable such behaviors. 

In the context of embeddability as a definition of dissipativity, we showed that one can always find 
behaviors that can be embedded in both a strictly dissipative behavior and a strictly 'anti-dissipative' 
behavior, thus raising a question on the embeddability definition. We related this question to the well- 
known result that the nullator one-port circuit is not realizable using only RLC elements. 

References 

[1] M. Athans, The role and use of the stochastic Linear-Quadratic-Gaussian problem in control systems 
design, IEEE Transactions on Automatic Control, vol. 16, pages 529-552, 1971. 

[2] V. Belevitch, Classical Network Theory, Oakland, CA: Holden-Day, 1968. 

[3] S. Bittanti, A. J. Laub and J. C. Willems (Eds.), The Riccati Equation, Springer- Verlag, 1991. 

[4] R. Bott and R. J. Dufiin, Impedance synthesis without transformers, Journal of Applied Physics, vol. 
20, page 816, 1949. 

[5] S. Boyd, L. E. Ghaoui, E. Feron and V. Balakrishnan, Linear Matrix Inequalities in System and 
Control Theory, SIAM, Philadelphia, 1994. 

[6] O. Brune, Synthesis of a finite two-terminal network whose driving-point impedance is a prescribed 
function of frequency. Journal of Mathematics and Physics, vol. 10, pages 191236, 1931. 

[7] M.K. Qamlibel, J. C. Willems and M. N. Belur, On the dissipativity of uncontrollable systems, in 
Proceedings of the 42nd IEEE Conference on Decision and Control, Hawaii, USA, December 2003. 

[8] H.J. Carlin, "Singular network elements," IEEE Transactions on Circuit Theory, vol. 11, no. 1, pages 
67-72, 1964. 

[9] L. E. Faibusovich, Matrix Riccati inequality: existence of solutions. Systems & Control Letters, vol. 
9, pages 59-64, January 1987. 

[10] I. Gohberg, P Lancaster, and L. Rodman, Indefinite Linear Algebra and Applications, Birkhauser, 
Basel, 2005. 

[11] A. G. J. MacFarlane, An eigenvector solution of the optimal linear regulator problem. Journal of 
Electronics and Control, vol. 14, no. 6, pages 643-654, June 1963. 

[12] D. Pal and M. N. Belur. Dissipativity of uncontrollable systems, storage functions, and Lyapunov 
functions, SIAM Journal on Control and Optimization, vol. 47, no. 6, pages 2930-2966, 2008. 

[13] H.K. Pillai and S. Shankar, A behavioural approach to control of distributed systems, SIAM Journal 
on Control and Optimization, vol. 37, pages 388408, 1998. 



25 



[14] J. E. Potter, Matrix quadratic solutions, SIAM Journal of Applied Mathematics, vol. 14, no. 3, pages 
496-501, 1966. 

[15] P. Rapisarda and J. C. Willcms, State maps for linear systems, SIAM Journal on Control and 
Optimization, vol. 35, no. 3, pages 1053-1091, 1997. 

[16] S. Karikalan and M.N. Belur, Uncontrollable dissipative dynamical systems, in Proceedings of the 
IFAC World Congress, Milan, Italy, 2011. 

[17] C.W. Scherer, The Riccati Inequality and State-space "Koa-optimal Control, Ph.D. thesis. University 
of Wiirzburg, Germany, 1990. 

[18] H. L. Trcntclman and J. C. Willems, Every storage function is a state function. Systems & Control 
Letters, vol. 32, pages 249-259, 1997 

[19] J. C. Willems, Dissipative dynamical systems - Part I: General theory. Part II: Linear systems with 
quadratic supply rates. Archive for Rational Mechanics and Analysis, vol. 45, pages 321351, 352393, 
1972, 

[20] J. C. Willems and H. L. Trentelman, On quadratic differential forms, SIAM Journal on Control and 
Optimization, vol. 36, no. 5, pages 1703-1749, 1998. 

[21] J.C. Willcms and H.L. Trentelman, Synthesis of dissipative systems using quadratic differential forms, 
IEEE Transactions on Automatic Control, vol. 47, pages 53-69, 2002. 

[22] J. G. Willems, Hidden variables in dissipative systems, in Proceedings of the 43rd IEEE Conference 
on Decision and Control, Bahamas, pages 358-363, 2004. 

[23] J. C. Willems, Dissipative dynamical systems, European Journal of Control, vol. 33, pages 134-151, 
2007. 

[24] J.G. Willems, The behavioral approach to open and interconnected systems. Control Systems Mag- 
azine, vol. 27, pages 46-99, 2007. 

[25] W. M. Wonham, Linear Multivariable Control: a Geometric Approach, Springer- Verlag, New York, 
1979. 



26 



